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ABSTRACT 

Aims. We present detailed modelling of the exoplanet-host star p Arae, using a new method for the asteroseismic analysis, and taking into 
account the new value recently derived for the Hipparcos parallax. The aim is to obtain precise parameters for this star and its internal structure, 
including constraints on the core overshooting. 

Methods. We computed new stellar models in a wider range than Bazot et al. (2005), with various chemical compositions ([Fe/H] and Y), 
with or without overshooting at the edge of the core. We computed their adiabatic oscillation frequencies and compared them to the seismic 
observations. For each set of chemical parameters, we kept the model which represented the best fit to the echelle diagram. Then, by comparing 
the effective temperatures, gravities and luminosities of these models with the spectroscopic error boxes, we were able to derive precise 
parameters for this star 

Results. First we find that all the models which correctly fit the echelle diagram have the same mass and radius, with an uncertainty of the 
order of one percent. Second, the final comparison with spectroscopic observations leads to the conclusion that besides its high metallicity, 
ji Arae has a high helium abundance of the order of Y = 0.3. Knowing this allows finding precise values for all the other parameters, mass, 
radius and age. 



1. Introduction 

Asteroseismology of solar type stars is a powerful tool which 
can lead to the precise determination of the stellar parameters 
when associated with spectroscopic observations. Special tools 
have been developped for this purpose, which help constraining 
the internal structure of the stars from the observable acoustic 
frequencies. In previous papers (e.g. Vauclair et al. 2008 ), we 
developped a systematic way of comparing models with obser- 
vations. First, for each set of chemical composition [Fe/H] and 
Y, we compute evolutionary tracks and for each mass the model 
which best fits the observed large separation is derived. Then 
we keep only the model that also best fits the observed echelle 
diagram. We finally have a set of models correctly fitting the 
seismic data, with different abundances. Interestingly enough, 
all these models have the same mass and radius, and of course 
the same log g (see Vauclair et al. l2008l for the star i Hor). All 
these models are placed in a log g - log Teff diagram, together 
with the observed spectroscopic boxes. Finally we only keep 
the models which satisfy all the constraints, including seismol- 
ogy and spectroscopy. In this way, precise values of the stellar 
parameters are found. 

Here we test this method on the exoplanet-host star fi Arae 
(HD 160691). This G3 IV-V type star is at the centre of a four- 
planets system: 
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- ij. Arae b, a 1.7 Mj planet with an eccentric orbit (e - 0.31) 
and a period P ^ 638 days (Butler et al. 2000 ^ Jones et al. 
I2OO2I 1. 

- /i Arae c, a 14 M© planet with P ^ 9.5 days. (Santos et al. 
,2004bJ . 

- yU Arae d, recently discovered by Pepe et al. ( I2007I I. with a 
period P ^ 310 days and a mass of 0.52 Mj. 

- yU Arae e, a long period planet (P ^ 1 1.3 years) with a mass 
of 1.81 Mj (Jones et al. l2002l McCarthy et al. l2004l Pepe 
etal. lIOOTl i. 

The star's overmetallicity has been established by many 
groups of observers: Bensby et al. (2003), Laws et al. (I2003I I. 
Santos et al. C2004 a), Santos et al. G2004bJ and Fischer & 
Valenti (l2005l l (see Table©. 

fi Arae was observed with the HARPS spectrometer at La 
Silla Observatory during eight nights in June 2004 to obtain ra- 
dial velocity time series. The analysis of these data led to the 
discovery of up to 43 frequencies that could be identified with 
p-modes of degrees ^ = to 3 (Bouchy et al. l2005l l. A detailed 
modelling was given by Bazot et al. (I2005I I. They computed 
models with two different assumptions that could explain the 
observed overmetallicity: overabundance of metals in the orig- 
inal interstellar cloud or accretion of planetary material onto 
the star. They tried to obtain evidence of the origin of yu Arae's 
overmetallicity. The results were not conclusive in that respect, 
as the differences were not large enough to decide between the 
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Table 1. Effective temperatures, gravities, and metal abun- 
dances observed for fj. Arae. [Fe/H] ratios are given in dex. 



Teff(^) 


logg 


[Fe/H] 


Reference 


5800±100 


4.30±0.10 


0.32±0.10 


Bensby et al. 2003 


58 11 ±45 


4.42±0.06 


0.28±0.03 


Laws et al. 2003 


5798±33 


4.31+0.08 


0.32+0.04 


Santos et al. 2004a 


5813±40 


4.25±0.07 


0.32+0.05 


Santos et al. 2004b 


5784±44 


4.30±0.06 


0.29+0.03 


Fischer &Valentil2005l 



two scenarii. Later on, other evidences were obtained that the 
observed overmetallicity in exoplanet host stars must be origi- 
nal, not due to accretion (see Castro et al. 2009 ). 

We computed new models, testing various values of the 
original metallicity and helium abundance. We compared the 
parameters of these models with those obtained from spec- 
troscopy, and introduced the new value of the Hipparcos paral- 
lax, as given by van Leeuwen (I2007I I. We also analysed models 
with overshooting at the edge of the stellar core. In some of 
these models, negative small separations appear so that there 
is a crossing point in the echelle diagram for the lines i = Q - 
{ -2. This effect was specially discussed in Soriano & Vauclair 
( I2OO8I 1. who showed that all solar type stars go through a stage 
where the small separations become negative in the observ- 
able range of frequencies. Here, we present a direct applica- 
tion of this theoretical effect, which is used to constrain core 
overshooting. 



2. Stellar parameters 

2.1. Spectroscopic constraints 



In Bazot et al. (I2005I I. the effective temperatures and luminosi- 
ties of the models were the first parameters used for compar- 
isons with the observations. We now prefer to compare the 
models with the spectroscopic observations in a more consis- 
tent way, by using the "triplets" Tgff, log g and [Fe/H]. These 
three parameters are indeed consistently given be the same ob- 
servers, while the luminosities are obtained in a different way, 
using Hipparcos parallaxes. Here we compare the luminosities 
as a second step. For this reason, we do not use exactly the 
same set of references as given by Bazot et al. j2005l) . as we 
only keep those which give constraints on the surface gravity 
(see Tableffll- 

2.2. Tiie luminosity of n Arae 

The visual magnitude of ji Arae is V = 5.15 (SIIVIBAD 
Astronomical data base). Van Leeuwen (I2007I I carried out a 
new analysis of the Hipparcos data. He derived a new value of 
the parallax of fi Arae: n - 64.48 ± 0.31 mas, which is lower 
than the first one derived by Ferryman et al. ( 1997) and has a 
reduced error. Using this new value, we deduced an absolute 
magnitude of Mv - 4.20 + 0.04. 

From Table [1] we have an average value of effective tem- 
perature for yu Arae of T^eff - 5800 + 100 K. Using the tables of 



Flower (fT996] l. we obtained BC = -0.08+0.02 for the bolomet- 
ric correction. With a solar absolute magnitude of Mhoi,& - 4.75 
(Lejeune et al. 11998 1). we deduced for ji Arae a luminosity of 
log (L/Lq)= 0.25 + 0.03. This value is lower than the one de- 
rived by Bazot et al j2005T l in their analysis of this star, namely 
log (L/Lo)= 0.28 + 0.0012. 

2.3. Seismic constraints 

The star ju Arae was observed with the HARPS spectrograph, 
dedicated to the search for exoplanets by the means of precise 
radial velocity measurements. HARPS is installed on the 3.6m- 
ESO telescope at la Silla Observatory, Chile. 

These measurements led to the identification of 43 p-modes 
of degrees = to ^ = 3 with frequencies between 1.3 and 
2.6 mHz (Bouchy et al. 120051 1. The frequency resolution of the 
time-series was 1 .56 yuHz, and the uncertainty on the oscillation 
modes has been evaluated to 0.78 //Hz. The mean large separa- 
tion between two modes of consecutive order, computed in the 
observed range of frequencies, is Avo - 90 //Hz, with an un- 
certainty of 1 . 1 //Hz. The small separations present variations 
which are clearly visible in the echelle diagram. Their average 
value in the observed range of frequencies is < 6v >= 5.7 //Hz. 

3. Evolutionary tracks and models 

We computed evolutionary tracks with the TGEC code (Hui 
Bon Hoa l2008l Richard et aL [T996] l. with the OPAL equations 
of state and opacities (Rogers & Navfonov 120021 Iglesias & 
Rogers 119961 ) and the NACRE nuclear reaction rates (Angulo 
et al. 119991) . For all models, the microscopic diffusion was in- 
cluded as in Paquette et al. (1986) and Michaud et al. j20041 ). 
The treatment of the convection is done in the framework of 
the mixing length theory and the mixing length parameter is 
adjusted as in the Sun: a - 1.2, (Richard et al. 120041) . We also 
computed cases with a - \.l and a - 1.9 to test the corre- 
sponding uncertainties. We found that for 1 . 10 M© and a - 1 .7, 
the track is very close to that of 1.12 M© and a - 1.8. On the 
other side, the 1.10 Mq and o- = 1.9 track is close to the 1.08 
Mq and a = 1.8 one. This was taken into account in the deter- 
mination of the uncertainties on the final results (Sect. 5). 

From the literature, we found two values for the metallic- 
ity of // Arae: [Fe/H] =0.29 (Laws et aL l2003l Fischer & Valenti 
I2D05 ) and [Fe/H] =0.32 (Bensby et al. l2003l Santos et al. l2004al 
and 2004b). For each value of [Fe/H], we first computed two 
series of evolutionary tracks for two different values of the he- 
lium abundance: 

- a helium abundance increasing with Z as given by the 
Isotov & Thuan j2004l ) law for the chemical evolution of 
the galaxies: we label this value Yq. 

- a solar helium abundance Y© = 0.2714 

The results obtained in the log g - log Tp// and log (L/Lq) - 
logTe// planes are displayed in FiglT] Each graph coiTesponds 
to one value of the metallicity. On each graph, the five error 
boxes are plotted, but those that coiTespond to the same value of 
[Fe/H] as given by the groups of observers are drawn in thicker 
lines. 
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We then computed two similar series including overshoot- 
ing at the edge of the convective core. Here overshooting is 
simply described as an extension of the convective core by a 
length a„y Hp, where Uov is the overshooting parameter, and Hp 
the pressure scale height. In these two series, the overshooting 
parameter is fixed at Ogv - 0.2 (Fig.|6]l. Finally, we tested more 
precisely the seismic constraints on overshooting by varying 
aov in small steps between 0.0 and 0.2 for models of masses 
1.1 Mo (Fig.|3. 



4. Seismic tests and results 

4.1. Models without overshooting 

We computed the adiabatic oscillation frequencies for a large 
number of models along each evolutionary track with the 
PULSE code (Brassard & Charpinet'2008'). These frequencies 
were computed for degrees = to ^ = 3, and for radial orders 
typically ranging from 4 to 100. For each track, we selected the 
model that has a mean large separation of 90 yuHz, computed 
in the same frequency range as the observed one. We plotted 
on each graph (Fig.[T]) the corresponding iso-Avo 90 //Hz line. 
The parameters of these models are given in Tables |2] to |5] We 
can check that all the models that fit the same large separation 
present the same value of the parameter g/R or M/R^ where g, 
M and R are respectively the gravity, mass and radius of the 
star. 

This is because the large separation approximately varies 
like c/R, where c represents the mean sound velocity in the star, 
which itself approximately varies like -s/tJJj. With the usual 
scaling for the stellar temperature T ^ i^M/R, we obtain that 
(c/Ry^ varies like M/R^, independently of the chemical com- 
position. 

Note on Fig. [T] that the iso-Avo 90 juHz lines never cross 
Laws et al. (I2003I I or Santos et al. (I2004al l error boxes. Models 
with Av - 90 yuHz have a value of log g much lower than the 
one derived by these two groups of observers. 

Kjeldsen et al (I2008I I discussed possible corrections on the 
frequencies due to near-surface effects. They gave a parame- 
terisation formula including three parameters, in which two of 
them can be deduced from the third one. As this parameter is 
unknown, they suggest using the same value for all solar type 
stars as deduced for the Sun in a first approximation. They find 
that the induced shifts increase with the frequency values, so 
that it can slightly modify the computed mean large separation, 
according to the frequency range. In our echelle diagrams for /i 
Arae, the models and the observations correctly fit when we do 
detailed mode to mode comparisons (Fig.|4|. We conclude that 
these near-surface effects are negligible in the observed range 
of frequencies, i.e. at least up to 2.3 mHz. 

For each value of [Fe/H] and Y we searched, among the 
models which fit a mean large separation of 90 //Hz, the one 
which also gives the best fit to the observational echelle dia- 
gram. In order to find this "best model", we performed a 
minimization for the echelle diagram of the models. Tables |2] 
to |5] show the parameters of some models close to the mini- 



mum ofx^- We used the expression 

. l|./ v.(/)-v..,(O f 
^ \ o-ots I 

where Vobsii) is the observational modulated frequency, v,/,(/) 
its theoretical counterpart, and cTobs the uncertainty of the ob- 
served frequencies, for which we used the value given by 
Bouchy et al. ( [2005l l: 0.78 //Hz. 

For all sets of chemical composition, the values present 
a minimum for models of mass M = 1.10 M©, radius 
R - 9.46x10'" cm, and gravity log g - 4.21 (Fig. |2]i. These 
models represent the four "best models", one for each couple 
([Fe/H], Y). Their echelle diagrams are displayed in Fig. [3] For 
a better visibility of the figures, we give the observed mode 
frequencies separately and the modelled ones as lines. It is re- 
markable that their masses, radii and thus gravities are identical 
and has important consequences, which will be discussed be- 
low. 

They all lie at the beginning of the subgiant branch. For Yg, 
for both values of the metallicity, the models computed along 
the evolutionay track have a convective core during the main 
sequence phase which disappears at the present stage, leaving 
a helium core with sharp edges. Models with a solar helium 
abundance Y© do not develop a convective core during their 
evolution, they have a helium core with sharp boundaries. 

For a fixed value of [Fe/H], the age of the best model de- 
creases when the helium abundance increases (compare respec- 
tively Tables|2]and|4]and Tables[3]and|5]l. This is due to the dif- 
ference in the evolutionary time scales according to the value 
of Y. The main sequence time scale increases for increasing he- 
lium abundance. As can be seen in Fig. [1] for a given value of 
[Fe/H] the effect of a smaller helium abundance is to move the 
models to lower effective temperatures. For a given value of Y, 
the models are cooler for a higher metallicity. 

We plotted these four models in a log g - log Tgff plane, 
together with the spectroscopic error boxes (Fig.|5]l. We clearly 
see that models with the lower helium abundance Yq (models 2 
and 4) lie outside the spectroscopic error boxes. Their effective 
temperatures are too low compared to the values derived by the 
spectroscopists. On the other hand, the models computed with 
a high helium content (models 1 and 3) lie right inside the error 
boxes: their external parameters correspond to those derived by 
spectroscopy. These are also the models with the lowest (see 
Tables |2] to |5]l. The luminosities of all these models agree with 
that derived from the Hipparcos parallax, namely log (L/L©) - 
0.25 + 0.03 (Sect. 2.2). The echelle diagram for the best model 
3 is also presented in Fig.|4]in more detail, with all the observed 
and modelled oscillation modes given separately. 

These models differ from the overmetallic model proposed 
by Bazot et al. 42005), which corresponded to a more massive 
star (1.18 Mq) on the main sequence. One of the reasons for 
this difference is related to the fact that in Bazot et al. (120051 ) the 
luminosity was the basic parameter used for comparisons, and 
that its value has been modified. The present analysis is more 
precise and consistent in the comparison with the observations. 

Another reason for the different results obtained here is re- 
lated to the helium content of the star. We find that only models 



4 M. Soriano and S. Vauclair: New seismic analysis of Arae 

Table 2, Characteristics of some overmetallic models with [Fe/H]=0.32 and Yc 



Mass 


Age (Gyr) 


^ogg 


log Teff 


log L/Lo 


R (xlO"' cm) 


M/R^ 


1 

X' 


1.08 


7.112 


4.2121 


3.7594 


0.2519 


9.40 


2.58 


1.700 


1.10 


6.318 


4.2149 


3.7644 


0.2770 


9.46 


2.58 


1.611 


1.12 


5.748 


4.2172 


3.7682 


0.2975 


9.52 


2.58 


1.668 


1.14 


5.387 


4.2200 


3.7667 


0.2967 


9.57 


2.58 


1.742 


1.16 


4.953 


4.2246 


3.7642 


0.2894 


9.61 


2.59 


1.915 



Table 3. Characteristics of some overmetallic models with [Fe/H]=0.29 and Yc 



Mass 


Age (Gyr) 




log Teff 


log L/Lo 


R (xlO"' cm) 


M/R^ 


1 

X' 


1.06 


7.916 


4.2064 


3.7571 


0.2371 


9.35 


2.58 


1.734 


1.08 


7.152 


4.2092 


3.7607 


0.2569 


9.41 


2.58 


1.650 


1.10 


6.367 


4.2152 


3.7653 


0.2801 


9.46 


2.58 


1.620 


1.12 


5.716 


4.2145 


3.7698 


0.3039 


9.52 


2.58 


1.737 


1.15 


5.177 


4.2215 


3.7672 


0.3006 


9.60 


2.58 


2.100 



Table 4. Characteristics of some overmetallic models with [Fe/H]— 0.32 and Y{ 



Mass Age (Gyr) log g log T^ff log L/Lq R (xlO^" cm) M/R^ 



1.08 


9.561 


4.2088 


3.7459 


0.1981 


9.41 


2.58 


1.663 


1.10 


8.678 


4.2120 


3.7491 


0.2152 


9.46 


2.58 


1.632 


1.14 


7.002 


4.2183 


3.7576 


0.2589 


9.56 


2.59 


1.670 


1.18 


6.073 


4.2230 


3.7562 


0.2636 


9.68 


2.58 


2.060 


1.20 


5.512 


4.2281 


3.7576 


0.2714 


9.70 


2.60 


3.572 



with a high helium value (Y=0.30) correctly reproduce all the 
observations. The helium abundance, which was kept as a free 
parameter by Bazot et al. (120051 1. is now well constrained. This 
determination of Y represents an important improvement. 

Interestingly enough the case of yU Arae is different from 
that of L Hor, discussed by Vauclair et al. (12008k The astero- 
seismic study of i Hor showed that its helium abundance was 
low (Y=0.255), comtrary to that of yu Arae. These two differ- 
ent results have interesting implications for the details of the 
nuclear processes which lead to the observed overmetallicity. 
This will be discussed in the conclusion. 

4.2. Constraints on core overshooting 
4.2.1. Models with aov=0.20 

The observational echelle diagram presents characteristic vari- 
ations in the small separations. We were particularly interested 
in the fact that around v = 2 mHz the lines f = and { - 2 
come close to each other, so that the small separation is smaller 
than average. In a previous paper (Soriano & Vauclair I2008I I 
we noticed that the small separations between { - and i — 2 
may vanish and become negative in stars at the end of the main 
sequence phase due to their helium cores. In this case there is 
a crossing point in the echelle diagram, with an inversion be- 
tween ( - and { - 2 modes at high frequencies. 



Here we wanted to check whether adding an overshooting 
layer at the edge of the convective core could lead to such an 
effect and reproduce the observed echelle diagram in a different 
way. We tried to find a model that could fit the observational 
echelle diagram and present a crossing point between the lines 
^ = 0- ^ = 2 at a frequency close to 2 mHz. 

We computed evolutionary tracks with an overshoot- 
ing parameter ffo,, - 0.20 for the two values of metallic- 
ity: [Fe/H] =0.29 and 0.32, the two values of helium abun- 
dance: Yg and Yq, and for masses ranging from 1.06 to 
1.20 Mq (Fig. |6]l. As for the cases without overshooting, we 
picked up the models with a large separation of 90 yuHz and 
draw curves of iso large separations. 

In all these models, a convective core develops, even for 
masses smaller than 1.10 M©. Models with masses ranging 
from 1.06 to 1.14 lie at the end of the main sequence. 
They have a high central helium abundance Yc, between 0.8 
and 0.9. Most of these models present a crossing point between 
the { - Q and 1-2 curves in the observed frequency range. 
However, as can be seen on the example shown in Fig.[8](upper- 
left panel), the computed lines do not reproduce the observed 
ones. 

Models with masses higher than 1.14 Mq correspond to 
younger stars, still on the main sequence. Their convective core 
is well developped, but their helium-content Yc is not high 
enough to induce a discontinuity in the sound speed profile 
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Mass 


Age (Gyr) 


log 4' 


log Tgff 


log L/Lo 


R (xlO'" cm) 


M/R' 


x' 


1.10 


8.559 


4.2114 


3.7530 


0.2319 


9.47 


2.58 


1.627 


1.12 


7.766 


4.2151 


3.7566 


0.2504 


9.51 


2.59 


1.670 


1.14 


6.734 


4.2179 


3.7607 


0.2717 


9.57 


2.58 


1.711 


1.16 


6.227 


4.2197 


3.7650 


0.2946 


9.63 


2.58 


1.762 


1.18 


5.802 


4.2250 


3.7590 


0.2727 


9.65 


2.60 


1.849 



and so to obtain negative small separations in the considered 
range of frequencies. There is no crossing point in their echelle 
diagram. These models do not fit the observational echelle dia- 
gram either, as can be seen in the example in Fig.[8](upper-right 
panel). 

These results show that the observed frequencies in Arae 
cannot be interpreted in terms of a crossing point between the 
{ - Q and 1-2 curves. This is consistent with the interpre- 
tation of the doublets observed in the t - 2 curve in terms of 
rotational splitting (Bouchy et al. '2005V In this framework, the 
upper part of this curve could not be misinterpreted as ^ = 
frequencies. 

As overshooting with ao,. = 0.20 is excluded, we decided 
to test smaller overshooting parameters to obtain a strong con- 
straint on the possibility of overshooting in this star 

4.2.2. Constraints on the overshoot parameter 

We computed new evolutionary tracks for models with 
1.10 M0, [Fe/H]=0.32, and Yc, gradually increasing the over- 
shoot pai-ameter: = 0.001, 0.002, 0.005, 0.01, 0.04, 0.10 
and 0.20. The results obtained in the log g - log T^/y plane are 
displayed in Fig |2l Here again we picked up the models with 
a mean large separation of 90 yuHz (iso-Av line on the tracks). 
We plotted their echelle diagrams and performed tests. The 
internal characteristics of four of these models are presented in 
Tableland their echelle diagrams are displayed in Fig.[8](mid- 
dle and lower panels). Models with aov - 0.001 and 0.002 are 
undistinguishable from models without overshooting. 

The more overshooting is added at the edge of the stellar 
core, the more the evolutionary time scales are increased, as 
can be seen in Fig.|2l The developpement of the convective core 
is increased during a longer main sequence phase. The models 
that have a mean large separation of 90 juHz are not in the same 
evolutionary stage, depending on the value of Uov, and they do 
not have the same internal structure. Models with aav< 0.002 
are at the beginning of the subgiant branch, as are models with- 
out overshooting. When increases, with a value between 
0.002 and 0.10, the models are in the phase of contraction of 
the convective core. And finally, for aov> 0.10, the models are 
at the end of the main sequence. This difference of evolutionary 
stages explains that the central helium abundance is lower for 
models with a higher overshooting parameter. 

We studied the evolution of the oscillation frequen- 
cies when aov increases. For small values of overshooting 
{aov - 0.001 and Uov - 0.002), there is no visible influence 



on the oscillations frequencies. For Uov - 0.005, the reduced 
is a little higher (x^ - 1 .66) than for the model without over- 
shooting. We can see on the echelle diagram (Fig. [8] middle-left 
panel) that the lines 1-^-1-2 are closer for high frequen- 
cies. For aov - 0.01, these lines cross at v = 2.6 mHz, as can 
be seen in the echelle diagram (Fig.[8j middle-right panel). The 

has increased: 1.73. For higher values of Uov {a„^> 0.01), 
the crossing point appears for lower frequencies decreasing for 
increasing aov (Fig- E] lower panels). Their value is high 
(2.45 and 3.04, respectively). These models do not fit the ob- 
servations. We conclude from these computations that the over- 
shooting parameter is small in // Arae, less than 0.01. 

Let us recall however that in these computations overshoot- 
ing is simply treated as an extension of the convective core. We 
have demonstrated that the seismic analysis of stars is able to 
give precise constraints on their central mixed zone. This does 
not exclude other kinds of mild macroscopic motions, plumes 
for instance, provided they do not lead to complete mixing. 

5. Summary and discussion 

This new analysis of the star /i Arae confirms that seismology 
is very powerful and can provide precise values of the stellar 
parameters with the help of spectroscopy to raise the final de- 
generacy between the best computed models. 

The procedure we use may be summarized as follows. 

- We need at least the identification all = 0, 1 and 2 obser- 
vational modes, and hopefully also = 3 as in the present 
paper 

- Detailed comparisons of observed and computed modes are 
necessary. Global values of the so-called large and small 
separations are not precise enough to derive these parame- 
ters. The comparisons of detailed frequencies are needed to 
take into account the fluctuations and modulations induced 
by the internal structure of the star, in particular by the cen- 
tral mixed zone. 

- In this framework, we first compare models and observa- 
tions in the log g - log T^.// plane. We may check that what- 
ever the chemical composition, models with the same large 
separation have the same value of the Mjt>? ratio, as pre- 
dicted by the asymptotic theory. 

- We then go further, by computing for various chemical 
compositions (Y, [Fe/H]) the model which fits the best the 
observed echelle diagram, using 2iX^ minimization process. 

- One of the most important results at this stage, which was 
already pointed out for the case of the star t Hor by Vauclair 
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Table 6. Characteristics of the overmetallic models of 1.10 Mq, [Fe/H]=0.32, for several values of the oovershooting parameter. 





Age (Gyr) 


log T,ff 


log (L/Lo) 




Yc 


r„./R* 


0.005 


6.599 


3.7565 


0.2484 


4.2162 


0.9482 


0.050 


0.010 


6.674 


3.7537 


0.2323 


4.2167 


0.9419 


0.055 


0.040 


6.838 


3.7499 


0.2177 


4.2160 


0.9243 


0.061 


0.100 


6.912 


3.7482 


0.2096 


4.2173 


0.8826 


0.066 



et al (I2OO8I 1. is that these best models, obtained for vari- 
ous chemical compositions, all have the same mass, radius, 
and thus gravity. With this method, mass and radius are ob- 
tained with one percent uncertainty. On the other hand, the 
age of the star still depends on the chemical composition, 
basically the helium value. 

- Then we compare the position of these best models with 
the spectroscopic error boxes in the log g - log T^jf dia- 
gram. This leads to the best choice of the chemical compo- 
sition. In particular, we now constrain the Y value, which 
represents an important improvement compared to previous 
studies. 

- Finally we check that the luminosity of the best of all mod- 
els is compatible with that derived from the apparent mag- 
nitude and the Hipparcos parallax. 

The parameters found for ju Arae are given in Table |2l The 
uncertainties have been tentatively evaluated by allowing a pos- 
sible increase of 0. 1 for each set of models and an uncer- 
tainty of 0. 1 on the mixing length parameter. They take into ac- 
count that all computations done with various sets of chemical 
parameters converge on the same model values. One must keep 
in mind however that these uncertainties do not include system- 
atic effects which would occur if stellar physics was strongly 
modified (new opacities, new nuclear reaction rates, new equa- 
tion of state, etc.). 

The new model is different from that given in Bazot et 
al. (I2005I I. Apart from the constraint on the Y value, the ba- 
sic reason is that in this previous paper the comparisons were 
only based on the stellar luminosity, which was misleading due 
to the Hipparcos parallax, which was later modified. Also the 
method used at that time was not as precise as the one we now 
use. Note that the scaling of parameters may lead to wrong re- 
sults for stars in which the seismic modes cannot be precisely 
identified . This is the case for example for the mass proposed 
by Kallinger et al. (I2008I I for the star fi Arae, which is much 
too lai-ge(1.23 Mq). 

We also performed an analysis of the size of the mixed core 
by testing the implications of overshooting on the mode fre- 
quencies. We found a strong constraint on the possibility of 
core overshooting, treated as an extension of convection: the 
size of this extension must be less than 0.5 % of the pressure 
scale height (overshooting parameter). This does not exclude 
other kinds of mild boundary eff'ects at the edge of the core, 
provided that they do not lead to strong mixing. 

At the present time, we were able to perform this deep seis- 
mic analysis on two solar type stars hosting planets, both with a 
large metallicity (about twice solar): l Hor (Vauclair et al l2008i l 



and yu Arae (this paper). In both cases, precise stellar parameters 
could be obtained. We found however an important diff'erence 
between these two overmetallic stars. 

In L Hor, the helium abundance is low, even lower than the 
solar value, in accordance with the helium value determined 
for the Hyades stellar cluster. As other observational parame- 
ters also coincide, we concluded that t Hor is an ejected mem- 
ber of the Hyades. The reason why the helium abundance is so 
low in these stars while the metallicity is high is still a mys- 
tery, although it certainly depends on the mass of the stars that 
polluted the original nebula. 

In yu Arae, on the other hand, the helium abundance is large, 
as expected from the usual laws for the chemical evolution of 
galaxies (Isotov & Thuan'2004). This star was formed in a neb- 
ula which suff'ered normal pollution from stars with propor- 
tional yields of helium and metals. 

Seismology can lead to precise values of helium abun- 
dances in solar type stars, where helium cannot be directly de- 
rived from spectroscopy. This represents a success, quite apart 
from all other results and constraints, and will be of importance 
for the study of the chemical evolution of our Galaxy. 
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Fig. 1. Evolutionary tracks in the log g - log Tg/f (left panels) and log (L/L©) - log T^/f planes (right panels), for the two values 
of metallicities found in the literature: [Fe/H]=0.29 and 0.32. The five error boxes shown are from Bensby et al. 2003 (asterisks), 
Laws et al. 120031 (crosses). Santos et al. l2004al (diamonds). Santos et al l2004bl (ti'iangles). and Fischer & Valenti 12005] (black 
triangles). The represented masses are 1.06, 1.08, 1.10, 1.12, 1.14, 1.16, 1.18, and 1.20 Mq. The horizontal thick line represents 
the iso-Avo 90 fiHz line, i.e. those models that have a mean large separation of 90 /iHz. 
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Fig. 2. Evolution of the )^ value with the mass of the models, between 1 Mq and 1 .20 Mq, for the four sets of chemical composi- 
tion: [Fe/H]=0.32, Xq (solid Une), [Fe/H]=0.29, Xq (dotted Une), [Fe/H]=0.32, Yq (dashed line), and [Fe/H]=0.29, Yq (dotted- 
dashed Une). 
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Fig. 3. Echelle diagrams for the best models with a metallicity of [Fe/H]=0.29 (upper panels) and [Fe/HJ=0.32 (lower panels), 
and an helium abundance Yc(left panels) and Y© (right panels). The lines represent the theoretical frequencies and the symbols 
their observational counterpart. Sohd hues and diamonds are for ^=0, dotted lines and triangles are for l=\, dashed hues and 
crosses are for t=2, and dotted-dashed Unes and asterisks are for ^=3. 
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Fig. 4. Echelle diagram for the best model with a metallicity of [Fe/H]=0.32 and a helium abundance Yq. Crosses and plus signs 
are for the observed frequencies, open circles and squares are for their theoretical counterpart. 
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Fig. 7. Evolutionary tracks for models of 1.10 Mq, with Yq, [Fe/H]=0.32, and a variable overshooting parameter The tracks 
are for aoy= 0, 0.001, 0.002, 0.005, 0.01, 0.04, 0.10, and 0.20. The horizontal thick hne represents the iso-Av 90/zHz hne. 
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Fig. 8. Echelle diagram for overmetallic models with overshooting: 1.06 Mq, [Fe/H]=0.32, Yc, and aov-0.20 (upper-left panel); 
1.20 Mq, [Fe/H]=0.29, Yq, and a'o,=0.20 (upper-right panel); 1.10 Mq, [Fe/H]=0.32, Yc, and a^.^O.OOS (middle-left panel), 
aoy=0.01 (middle-right panel), a'o,,=0.04 (lower-left panel), and aoi ^O.lO (lower-right panel) . The symbols are the same as in 
Fig.E] 



